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Electronic Thermal Conductivity Measurements in Intrinsic Graphene 
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The electronic thermal conductivity of graphene is of fundamental interest and can play an important role 
in the performance of graphene nano-scale devices. We measure the electronic thermal conductivity, K e , in 
suspended graphene in the nearly intrinsic regime over a temperature range of 20 to 300 K. We present a 
method to measure K e using two-point DC electron transport at low bias voltages, where the electron and 
lattice temperatures are decoupled. We find K e ranging from 0.5 to 11 W/m.K over the studied temperature 
range. The data is consistent with a simple model in which heat is carried by quasiparticles with the same 
mean free-path and velocity as graphene's charge carriers. 



The thermal properties of graphene are under intense 
investigation^— , and they are of fundamental interest to 
access information about the interactions of phonons and 
electrons in graphene^—. The large heat conductivity of 
graphene offers exciting prospects for heat management 
in microelectronic devices 7 -. While there have been sev- 
eral experimental reports of the phononic thermal con- 
ductivity, K p , in graphene^ 7 -, there is no report of the 
electronic heat conductivity, K e , in suspended graphene. 
This is because in most regimes K p is much larger than 
K ei which makes it difficult to measure the amount of 
heat carried by the charged quasiparticles (electron and 
holes). Measurements of K e are needed to complete 
our understanding of the various electron-phonon cou- 
plings in graphene ^ 12 ' 16 , and graphene's electro-optical 
propertie d 5 ^ 17 . Moreover, measuring and controlling K e 
could have applications in heat management in heavily- 
doped nm-scale devices where K e could be comparable 
to or larger than Kp^. 

We present a carefully calibrated method to extract K e 
in graphene using DC electron transport in suspended de- 
vices. The accuracy of the method is dependent on high- 
mobility (annealed) devices. We measured five samples 
which showed consistent data, and focus our presentation 
on the two samples for which we have a complete data set. 
The extracted K e are compared with calculated values, 
K e - t h, for a diffusing gas of Dirac quasiparticles. The 
agreement between theory and measurements is quanti- 
tative for one sample and within a factor of two for the 
second device over the entire temperature range (20 - 300 
K) studied. Throughout the text we use T to designate 
the lattice (cryostat) temperature, and T e for the electron 
temperature. At very low bias, \Vb\ < 1 mV, T = T e . 
We first describe our devices, secondly we present our 
thermometry method to monitor T e , then show how we 
apply a controlled ATs using Joule heating, and finally 
extract K e from transport data. 

Figure 1 (a)-(b) shows tilted SEM images of Sample 
A and B respectively. Sample A is 650 nm long, 675 
nm wide, and suspended 140 ± 10 nm above the sub- 
strate (AFM measurement) which consists of 100 ± 2 
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FIG. 1. (Color online.) Suspended graphene devices, (a) and 
(b) tilted SEM images of 650 nm and 400 nm long suspended 
graphene transistors (Samples A and B). (c) and (d) R of 
Samples A and B versus gate voltage, Va, at T e = T = 11, 
50, 100, 150, 210 and 300 K, and V B = 0.5 mV. 



nm of Si02 (ellipsometry measurement) on degenerately- 
doped Si which is used as a back-gate electrode. Sam- 
ple B is 400 nm long, 1.05 /jm wide, and suspended 
175 ± 10 nm above a 74 ± 2 nm Si0 2 film on Si. To 
prepare the samples, we used exfoliated graphene and 
standard e-beam lithography with evaporated contacts 
Ti(5nm) /Au(80nm). The samples were suspended with a 
wet BOE etch such that their only thermal connection is 
to the gold contacts. We annealed the devices using Joule 
heating in situ by flowing a large current in the devices^ 
(up to 540 and 840 /LtA for A and B) while under high 
vacuum, 10 -6 Torr. Panels (c) and (d) show two-point 
resistance data, R, for Sample A and B respectively after 
annealing, versus gate voltage, Vg, which controls charge 
density. From the width of the R maximum at 11 K, we 
extract a half-width-half-maximum, HWHM, of 0.45 and 
0.6 V for Samples A and B. Using a parallel plate model 
for the gate capacitance of the devices, these HWHMs 
correspond to an impurity induced charge density 2 ^ of 
n* « 1.5 and 1.7 xlO 10 cm" 2 . 
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FIG. 2. (Color online.) Electron thermometry, (a) Temper- 
ature dependence of the two-point R in Sample A (black cir- 
cles, left axis) and Sample B (red squares, right axis) near the 
charge degeneracy no = 5.7 and 2.9 xlO 9 cm -2 . The dashed 
lines are numerically interpolated curves used for thermome- 
try. Inset: I — Vb data for Sample A, \Vb\ < 1 mV, whose 
slope is used to extract R. The solid (dashed) lines is at 
T = T e = 11 (300) K. (b) Relative conductance G/Ghk of 
Sample A (black circles) and Sample B (red squares) versus 
T = T e . The filled symbols show the raw two-point data, and 
the open symbols the data after subtracting R c = R -Dirac 
(see text). The solid and dashed lines are power law fits con- 
sistent with charge impurity scattering, (c) Relative conduc- 
tance in the hole doped regime, na = — 1.77 and - 2.0 xlO 11 
cm -2 , with the same symbols as in (b). Sample A shows 
diffusive transport, while Sample B shows ballistic transport. 



The devices were fabricated with large contact areas 
between the gold electrodes and graphene crystals, 1.2 
to 3 /im 2 per contact, to minimize the total contact re- 
sistance, R c . An upper bound for R c can be extracted 
from the two-point R — Vg curve in Fig. 1 (c) by fitting 
the data with the expression 21 R = R + (L/W)(l/nGe/J,) 
where R is the resistance due to neutral scatterers plus 
R c , L is the length of the device, W the width, uq the 
charge density induced by Vg, \i the mobility, and e the 
electron's charge. We fit the data at T = T e = 100 K 
for (Vg — Vb) > 1.3 V to avoid the thermal smearing 



around the Dirac point, Vb. The extracted mobility for 
Sample A in the doped regime is /z » 8.5 x 10 4 cm 2 /V.s 
at 100 K, and R a w 682 ± 1 and 1135 ± 4 Q for hole 
and electron doping respectively. The difference between 
hole, R -h and electron doping, i? _ e , is understood as 
an additional p — n barrier for the electron due to p- 
doping from the gold electrodes^ 1 -. At the Dirac point, 
we let Ro-Dirac = [Ro-h + Ro- e )/2 = 908.5 n. We 
note that R -Dirac is much smaller than the total re- 
sistance of Sample A, therefore R c < R - Dirac has at 
most a modest impact on our measurements in this de- 
vice. Sample B is shorter than A and as we will show 
below, it enters the ballistic regime when doped. It is 
therefore not possible to extract its R . The contact ar- 
eas of Sample B are larger, and its width wider, than 
for Sample A. Assuming a similar resistance per unit 
area as for A, we expect R c < 657 fl for B. Based on 
the reported thermal conductance of Au/Ti/ Graphene 
and Graphene/Si02 interfaces^, the thermal resistance 
of our contacts are several orders of magnitude lower than 
the one we measure below for graphene. Thus, the ther- 
mal resistance of the contacts can safely be neglected. 

Figure 2(a) shows R vs cryostat temperature, T, cali- 
bration curves for Sample A (black circles, left axis) and 
Sample B (red squares, right axis) near Vg = Vb- R 
data are extracted from the slope of the I — Vb data, as 
shown in the inset of panel (a) at 11 K (solid) and 300 K 
(dashed), for ± 1 mV bias where no Joule heating effect 
is present (T e = T). The data are taken at Vg = 0.5 V 
close to Vb = 0.33 V for Sample A, and at Vg = V 
(Vb = -0-1 V) for Sample B, corresponding to ug = 5.7 
and 2.9 x 10 9 cm -2 . The T dependence of the data shows 
an insulating behavior up to 200 K for Sample A and up 
to 300 K for Sample B. The dashed lines in panel (a) are 
interpolated and will be used as thermometry curves to 
monitor T e . Note that the thermometry is most accurate 
where the curves are steepest. 

Figure 2(b) shows the relative conductance 
G(T) / G\\k in the intrinsic regime extracted from 
panel (a). The T dependence of G in graphene, at low 
charge density, is strongly dependent on the type of 
charge transport. For ballistic transport, we expect a 
very weak temperature dependence at low T, and a linear 
dependence when ksT >> Ef 23 ■ In the diffusive regime, 
the expected temperature dependence depends on the 
type of charge scatterers, and G(T)/G\\k oc T a with 
a = -1, 0, 2 for acoustic phonon, short-range (neutral), 
and long-range (charged) scatterers respectivel y 16 ' 24 . 
The temperature dependence of real samples is expected 
to combine all three types of scattering. We fit the 
data with a function G/Ghk = 1 + AT P , and extract 
p = 1.85, 1.74, 1.72 and 1.63 ± 0.03 for Sample A with 
R c = Ro-Dirac and (open and filled circles), and Sam- 
ple B with R c = 657 and (open and filled squares). 
This T-dependence strongly supports diffusive charge 
transport dominated by long-range charge impurities, 
as reported in previous experiments on high-mobility 
device o 16 ' 19 and expected theoretically^. The small 
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FIG. 3. (Color online.) Electron heating, (a) R vs Vb at 
various T for Device A at Va = 0.5 V s=s Vd- (b) Zoom-in of 
the low Vb range for the 100 K data. Joule heating due to Vb 
raises the flake's average T e above T. T e is extracted using 
Fig. 2(a). Panels (c) and (d) show T e vs Vb in Sample A over 
a broad range of Vb, and at low Vb respectively. 

departure from a T 2 dependence is expected as the 
samples are not exactly at the Dirac point. 

Panel (c) shows the contrasting T dependence of con- 
ductance in the hole doped regime for Sample A (cir- 
cles) and B (squares). The expected temperature depen- 
dence in the metallic regime for long-range scattering is 
G/Ghk = 1 ~ BT 2 , and we find an exponent of 2.13 ± 
0.18 and 1.93 ± 0.16 for Device A with the two limiting 
cases of R c . The behavior of Device B is very different, 
with an increasing G versus T, indicating that the sample 
is in the ballistic regime at this doping 2 ^. We note that 
the T-dependence in Sample A approaches T _1 above 
200 K, in the metallic regime, as the electron-phonon 
scattering becomes more important. We conclude that 
both samples are in the diffusive regime at low charge 
density, Fig. 2(b) and Supplemental Material (SM) sec- 
tion 1, and scattering is predominantly due to charged 
impurities. The well behaved data in Fig. 2(a), and its 
agreement with theory, serves as a reliable thermometer 
for T e in our devices. 

After establishing the thermometry for T e , we demon- 
strate controlled Joule self-heating of the electrons to ap- 
ply a thermal bias AT = T e — T between the suspended 
graphene and the gold electrodes. Figure 3(a) shows R 
vs Vb for Sample A over a broad range of Vb (data for 
Sample B is shown in Fig. SI). Panel (b) shows the 
details of data at Vb < 30 mV. R decreases monotoni- 
cally with increasing Vb, at all T. Using the calibration 
curve, Fig. 2(a), for the samples, it is possible to de- 
termine the average T e in the suspended device vs Vb, 
as shown for Sample A in Fig. 3(c) and (d). At very 
small Vb, T e increases superlinearly, then linearly be- 



tween ss 20 and 100 mV, and then sublinearly for Vb > 
100 meV. The transition from linear to sublinear T e vs 
Vb was previously observed in noise thermometry of T e in 
graphene^, and self-heating control of T e at low Vb was 
also demonstrated 2 ^. We restrict our K e measurement to 
Vb ^S, 30 meV, where the thermometry is most reliable. 
This rules out any R change due to scattering from opti- 
cal in-plane phonons, w 200 meV, and flexural phonons, 
« 70 meV, in graphene 16 . Phonons in the substrate can 
also be ruled out as the samples are suspended. Fig. 3(d) 
shows that the accuracy with which T e can be measured 
is much better than 1 K. We estimate T e errors from the 
scatter of the data in panel (d) , and similar plots at each 
T, to vary from 0.3 K (steepest regions of Fig. 2(a)) up to 

2 K (fiat regions of Fig. 2(a)). The smooth dependence 
of T e on Vb at all T is consistent with electrons having 
a well defined temperature, as predicted by calculations 
of the e — e collision length 2 -. This is also confirmed by 
the K e data shown below. 

Since our devices are wide, the effect of their edges on 
transport should be small. We use a 1-d heat equation 
to extract K e in our devices; K e ^^f + Q = where Q = 
RI 2 /WLh is the Joule heating power per unit volume, 
W the width of the flake, L the length, and h = 0.34 
nm the thickness. Using boundary conditions T e = T 
at the two ends (contacts) of the flake, we find T e (x) — 
T+(LQx— Qx 2 )/2K e . Averaging over the length we find, 

T e = (1/L) J^T e (x) dx = T+ (QL 2 )/{12K e ). Finally, 

K e = j§^, where AT = T e -T. Using R and I from Fig. 

3 and similar plots, for AT = 1, 2 and 5 K, we calculate 
K e vs T e as shown in Fig. 4 (a) and (b) for Samples A and 
B. Data in Fig. 4 (a)-(b) show a strong K e dependence 
on T e ranging from roughly 0.5 W/K.m at 20 K to 11 
W/K.m at 300 K. The T e range is limited to the region 
where we have accurate thermometry (Fig. 2(a)), up to 
200 K for A and 300 K for B. Error bars representing the 
accuracy of our thermometry 2 ^, are shown only for the 
AT = 5 K data for clarity. If the Vb needed to apply 
AT dopes significantly the sample, it could affect the 
measured K e . Using ntot{T) (SM section 1), we define an 
effective chemical potential /i e //(T) = hvFy/irntotiT). 
For instance, at T = 100 K, /i e // = 18 and 18.5 meV 
respectively for the two samples. The Vb necessary to 
achieve AT = 2 K at 100 K in Fig. 4 for Samples A 
and B are 12.1 and 10.0 mV, significantly smaller than 
Heff- We only observe a change in the extracted K e 
values when AT exceeds 8 K, and Vb > faff- Thus Vb 
does not affect our measured K e , with the caveat that 
we cannot measure K e precisely at n = 0. 

We compare our data with the usual model for diffus- 
ing particles in 2-dimensions, K e - t h = \Cvl. If the heat 
flow is due to charge carriers, then the specific heat is 
C = C e , the velocity is vp = 10 6 m/s, and the mean 
free path I is the same as for charge transport. We find 
(SM section 1), lA-avg and Is-avg = 71 ± 3 (85) and 47 
± 7 (59) nm on average over the T e range studied with 
R c = {Ro-Dirac)- We calculate C e using the density 
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FIG. 4. (Color online.) Electronic thermal conductivity, K e , 
in the quasi-intrinsic regime, n to t,T=o ~ 1.7 xlO 10 cm" 2 , (a) 
and (b) K e vs T e for AT = T e - T = 1 K (blue triangles), 2 K 
(black circles), and 5K (red squares) for Samples A and B. The 
solid lines are calculations of K e - The dashed lines show the 
same calculations with a contact resistance R c = R -Dirac- 
The error bars are only shown for the AT = 5 K data for 
clarity, and represent the errors in thermometry. 



of states for graphene and the Fermi-Dirac distribution 
(SM section 3). We plot K e -th in Fig. 4(a)-(b) as solid 
lines, R c — 0. They capture the qualitative behavior 
of our K e data, and are in good quantitative agreement 
for Sample A. The dashed lines shows K e -th if we use 
Rc = Ro-Dirac- If we account for R c , i.e smaller Q, K e 
changes by the same magnitude as K e _ t h but in the op- 
posite direction (not shown for clarity). The quantitative 
agreement between data and theory is not as definite for 
Sample B since its R is smaller than for Sample A, and 
the impact of R c may be important. The data and cal- 
culations shown in panel (b) with R c — are within 30 
% of each other, and if we include the maximum R c = 
657 f2, they are within a factor of two while preserving 
the same T e dependence. We fit a power law expression 
K e cx T p over T e = 45 - 185 K for both samples, and find 
p = 1.73 ± 0.15 and 1.63 ± 0.13 which is very close to 
the fit on K e -th, Pth = 1-62 and 1.59. As expected pth 
goes to 2 when (j,/kT <C 1. We conclude that the mea- 
sured thermal conductivity is consistent with the heat 
being carried by particles moving with the vf and I of 
the charge carriers. The magnitude of K e reaches « 11 
W / K.m at 300 K and n to t,T=o ~ 1-7 xlO 10 cm" 2 . 

A condition to make reliable K e measurements is that 
all of the Joule heat remains in the carriers until they dif- 



fuse to the leads. Both experiments and theory confirm 
that the electron-phonon energy transfer in high mobility 
graphene, at low Vb, is very small below 300 K£ ' 17 i 24 i 25 , 
and decreases at lower T and n. In our devices, we ex- 
tract a cooling length for hot electron (SM section 4), £ ss 
100 to 10 urn for T e = 20 to 300 K. Since £ is much larger 
than the length of our devices, and Vb below the energy 
of optical phonons, we expect T e and T to be decoupled 
in our devices when Vg ^ 0, and all of the Joule heat 
to be carried to the contacts by charge carriers. Indeed, 
the K e we measure are two to three orders of magnitude 
lower than the reported phonon thermal conductivity K, p 
in graphene^ 7 -. 

In summary, we fabricated high quality suspended 
graphene devices, developed self-thermometry and self- 
heating methods to measure and control T e , and made 
detailed measurements of the electronic thermal conduc- 
tivity in graphene. We measured K e in the quasi-intrinsic 
doping regime, n to t,T=o ~ 1-7 xlO 10 cm" 2 , from T e = 20 
K to 300 K. The measured values are in good agreement 
with a simple model of diffusing Dirac quasiparticles. 
This work was supported by NSERC, CFI, FQRNT, and 
Concordia University. We made use of the QNI clean- 
rooms. We thank Andrew McRae for discussions. 
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